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O ' Abstract 
(N 

\ We calculate the Wigner distribution function for the Calogero-Sutherland system which 

consists of harmonic and inverse-square interactions. The Wigner distribution function is sep- 



in 



arated out into two parts corresponding to the relative and center-of-mass motions. A general 
expression for the relative Wigner function is obtained in terms of the Laguerre polynomials by 
introducing a new identity between Hermite and Laguerre polynomials. 
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^ . PACS: 03.65.-w Quantum mechanics, 03.65.Sq Semiclassical theories and applications, 05.30.-d 

' Quantum statistical mechanics. 
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t3 ' 1 Introduction 
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O" 1 . Since the introduction of the Wigner function (WF) in 1932 for inclusion of quantum corrections 
>- ! to classical result s[l], phase space representations of the quantum mechanics have been a focus 
of continuing interest and found a wide range of applications [2 j. In this formalism one defines a 
distribution function W (q,p) of position q and momentum p in such a way that to every normalized 
state vector if> there corresponds a distribution function. In order to define the same physical system, 
W should be a Hermitian form of ip , and if W is integrated over p it should give the proper 
probabilities of the different values of q or vice versa. Since the WF is a probability distribution 
function one expects, as a natural condition on W(q,p), that it should be non-negative for all values 
of q and p: W (q,p) > 0. However , Wigner [3] proved that this result is incompatible with the first 
two conditions and it is now a common practice to work with a distribution function taking negative 
values for certain q and p in the phase space. 

The quantum Calogero-Sutherland model (CSM) having a quadratic confining q 2 plus an inversely 
quadratic 1/q 2 potentials|4j has applications in a wide variety of different areas of many body physics, 
due to the connection of its variants and itself directly with the hierarchical fractional quantum Hall 
effect jS], free oscillators on a circle [S], the spectrum of the Chern-Simons matrix model [7], short range 
Dyson model[8], and Witten-Dijkgraff-Verlinde equation[9j. Additionally, many works have also been 
realized to construct its N-fermion version [TU], algebra unification [TTJ, shape invariance[12j. gen- 
eralized statistics [13], statistical properties of quantum quasi-degenaracy[H], equivalence to decou- 
pled oscillators [T5]. 
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In addition to allowing one to analyze the dynamics of quantum systems entirely in phase space 
and thus to make comparison between their classical and quantum evolutions, there is also experi- 
mental interest on the measurements of WFs for certain quantum systems to probe the predictions 
of quantum mechanics, since WF contains complete quantum mechanical information as the wave 
function or density matrix has. A number of experiments have been reported where the measure- 
ment of WFs carried out for both vacuum and quadrature-squeezed states of light [16], molecular 
vibrational states [17]. various quantum states of the motion of a harmonically trapped atom [18], as 
well as for a massive particle wave packet [19J. Of particular importance are the experiments upon 
which the negatives in WFs corresponding to Fock states[H] and a superposition of macroscopically 
separated parts of matter field [19] have been observed. In this regard, due to the fact that it is now 
possible to realize new quantum mesoscopic devices such as quantum dots and quantum antidots 
by various experimental techniques, the CSM may serve as a model of two non-interacting electrons 
with an individual quantum antidot confined in a quantum wire or a stripe, where the repulsive 
inverse-square quantum antidot potential acts as a scattering center for electrons, or may be used 
as a one dimensional exactly soluble band model or quantum dot arrays wherein the coupling con- 
stant of the inverse-square potential of the CSM is chosen as > g > — 1/2 in dimensionless units, 
i.e., attractive [20l [21]. Moreover, very recently, Li et al. [22] have solved CSM with pseudo-angular 
momentum operator method, and they showed that, by discussing its several variants, the radical 
equations of three dimensional isotropic oscillator and hydrogen-like atom in both spherical and 
parabolic coordinates, one dimensional three body problem and the s-state of Morse potential all 
reduces to CSM. Therefore, WFs of CSM may provide a solid basis for the discussions of transport 
properties [23] of the above mentioned nanostructures. With these motivations, we study the WFs 
of CSM which has not only a particular significance in itself, but also enables us to understand the 
phase space picture of its variants and itself as well. 

The WF for two particles is defined by 



i r+co r+co 

W(qx,q 2 ;pi,P2) = — -5 / dyi / dy 2 

(7m) J-oo J~oo 

X** {q x + y x , q 2 + y 2 ) * (<?i - Vi, q 2 ~ Vv) 

x exp [2i {pxyi + p 2 y 2 ) /h\. (1) 

If we express the WF in terms of the center-of-mass and relative coordinates through the relations 

qi + q 2 = 2Q , qi -q 2 = q 
pi+p 2 = IP , Pl -p 2 = p 
2/1 + 2/2 = 2Y , yi -y 2 =y 

then, provided that the interparticle potential depends on the relative coordinates, equation ([I]) 
becomes 

i /*+oo /»+oo 

W(q,Q;p 1 P) = — § dy / dY 

(7m J J-oo J-oo 

xV*(q + y) (Q + Y)^{q- y) *(Q - Y) 

x exp [i (APY + py) /%) , (2) 

since the solution of the corresponding Schrodinger equation can be represented by a product of two 
functions, one for the center-of-mass and other for the relative coordinates. Then, Eq. ([2]) can be 
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separated as and relative WFs, W(q, p) and W (Q, P), respectively. We can therefore define the WFs 
for the center-of-mass and relative motions in the form 



i p+oo 

W(Q,P) = — dY ^* (Q + Y) *(Q - Y) exp [4iPY/fc] , (3) 

1 p+oo 

W{q,p) = — dy ^* (q + y) * (q - y) exp [ipy/fc] , (4) 

respectively. It should be noted that, while \l/'s in Eq. ([1]) represent the two body wavef unctions, 
\l/'s in Eq. (T5]-4) are single particle wave functions. Now, it is possible to present a coupled system of 
linear partial differential equations corresponding to the above defined WF, which requires the direct 
computation without solving the wave functions |24j. The Wigner representation is very convenient 
for studying quantum systems with Hamiltonians that include quadratic coordinates and momenta, 
since in this case the Wigner distribution function represents a good approximate description of the 
dynamics involved. The method, however, is not easy to handle when the potential contains higher 
order powers of coordinates, since this case comprises a differential equation for the Wigner function 
with terms as much as the number of the order . In recent years there has been a number of works to 
calculate the WF for various type of potentials: Infinite square well [25]. a double well potential [26J, 
the Posch- Teller potential[27j, the Morse oscillator [28]. a quantum damped oscillator [29|, the hy- 
drogen atom|30J, the rotational motion of a spherical top [31] are notable applications. A discrete 
WF for non-relativistic quantum systems with one degree of freedom has been developed in finite 
dimensional phase space and applied to a few simple system[32j. 

The layout of this paper is as follows: In Sec. [21 we discuss WFs for the center-of-mass and 
relative motions, and obtain general expressions in terms of Laguerre polynomials. In Sec. [3X we 
derive a new identity between Hermite and Laguerre polynomials to obtain a compact form for the 
WF of the relative part, and plot some of them for a few states to give an idea on their phase space 
behaviors. 
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The Hamiltonian describing two particles interacting pairwise by the Calogero- Sutherland potential 
is given by 

h 2 d 2 

=- 4 

2 



2m dq\ 



1 1 2 

+ mwlq 2 i +-^U{\q i -q j \) 



(5) 



where the second term is the confining potential and 

U - qj\) = [muj 2 (qi - qrf + 2g / (q { - qj) 2 ] /2 

simulates further interactions between two particles. If we now use the above defined center-of-mass 
and relative coordinates, then the relevant Schrodinger equation is separated out as a center-of-mass 
equation, which is a ID harmonic oscillator equation 



ft" d 2 1 

'2M^ + 2 M ^ 



*(Q) = E Q *(Q) 



(6) 
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and the Calogero-Sutherland system 



d 2 1 2 2 , # 



v& (g) = £7 9 * (q) , (7) 



where we have defined u 2 = oj 2 + 2u;q as a hybrid frequency. In Eqs. (jHl)-©, M and \i are total and 
reduced masses, and are given by 2m and m/2, respectively. 

The Wigner function corresponding to the center-of-mass motion defined by Eq. ([3]) through the 
solution of Eq. (jSJ) is well-known, and given by [2] 

W t (Q, P) = ^-J- exp \-Muj.Q 2 lh - P 2 /Muj.h 

V / 7m L 

, ,'2Mu;._ 2 2P 2 

xL e — — Q 2 + 



h hMuj, 

where P = 2P is used for the sake of comparison with the results presented in the associated 
literature, and £ takes values 0, 1,2,.. .. Although it is possible to obtain the WFs for the simple 
harmonic oscillator (Eq. ([8])) in various ways, for example, by using algebraic methods or by solving 
ordinary differential equations of WF |33j, the WF corresponding to the relative motion resulting 
in the CS system cannot be obtained by either methods. Therefore, we are compelled to obtain 
the corresponding Wigner function through solving Eq. (J7J). First, we make a change of variable by 
z = (iiuj /h) l l 2 q , which transforms Eq. ([7D into 

+ Un + 2/3 + 2 - z 2 + 1/4 * = 0, (9) 

where the new parameters are given by 

E q = hjj (2n + (5 + 1) , 1/4 - P 2 = -2fig/h 2 

with n = 0,1,2,... and g > —h 2 /8fi. It should be noted that these are the energy levels of 
one-dimensional isotropic harmonic oscillator with odd quantum numbers shifted by an amount 
(/3 — 1/2) hut. The solution to Eq. (Q can then be written [31] in terms of the Laguerre polynomials 

*n (?) = C n b a / 2 q a exp [-bq 2 /2] L^ 2 (bq 2 ) , (10) 

where a = (3+1/2, b = fiu/h and the normalization constant is given by C n = b 1 / 4 [n\/T (n + a + 1/2)] 1/2 . 
We can now build the associated WFs for the relative motion (RM) with the wave functions given 
by Eq. (1101) according to the definition of Eq. (J3J), which results in 

\(~^ \ 2 r+oo 

W na (q,p) = l -^J-b«exp[-bq 2 ] J dy (q 2 - y 2 )" 
x exp [-by 2 ] LT* [b (q + yf] 

xLn~* [b(q - yf] exp(2ipy/H), (11) 
where p = 2 p is used. If we use the binomial expansion of 



a + d) a = J2( ( * \a a -P d?, 

13=0 V ^ / 
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and series expansion of the Laguerre polynomials 



L k , x \ = f k + n \ ,„ 

n Z-^i m \ \ n — m J 

m=0 v 7 



X 



then Eq. (ITT1) is expressed in the form 



W na (q,p) = ^b a exp l-bq 2 ] 



m!r! \ n — m j \ n — r 

m r 



x £ ("If ( «0 ) g 2 - 2 ^ ( 2 ™ ) g 2 — 

/3=0 p=0 V ^ / 

2r 

X 



£(-!)'( 2 ; ) g 2 ^" F 2/3+M+P (p) , (12) 



where F 2j g +/ , +p (p) is given by the following integral 

/+oo 
cfy y 2 ^? exp [-% 2 + 2ipj//ft] . (13) 
-oo 

It is easy to show that this last integral can be expressed in terms of the Hermite polynomial as 
follows: 

F n (p) = ^ (n+1)/2 F ^ f exp [-pW] H n ( JL) , (14) 

with n = 2/3 + p + p[35j. Furthermore, the use of the relation[35j H n (u) = (— l) n e" 2 <9™ e - " 2 for the 
Hermite polynomials reduces the above expression to 

|(-f (15) 

where <9p represents the differentiation with respect to p. Hence, the WFs for the RM becomes, with 
these new definitions, 

W 7 ™ (g,p) = —j=- exp [-6g 2 ] 

V 7T0ft 



E^F( a ;ir)-( 9 -f%) : 



/'-•n™ / i , „ \ / \ 2m 

x 

m=0 

^ (-l) n /a-l + n V r / ih ^ 

r=0 v 7 v 



(16) 
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By using once again series expansion of the Laguerre polynomials, it is possible to express Eq. ffT6]) 
in an implicit form as well 



W na (q, p) 



nl 



y/^hT (n + a + l/2) 



:LT 1/2 



L ^ {g - ih d p ) 2 



x 



^ („2 I k2r& 

L h 



(a 2 + h 2 dl) 



ra-l/2 
p 2 /ifiwh 



— {q + ih d p ) 



(17) 



If we now define a unit of dimension by I = ^/h/muj, with /i = m/2 and M = 2m, then we can make 
positions and momenta dimensionless by Q = Q/l , P = IP/h, q = q/l , p — Ip/H , and frequency 
by uJ = uj/uj,. Hence, the relevant WFs for the center-of-mass and relative motions, Eqs. (jHJ) and 
( 1T71) . become 



W e (Q,P) 



-l) e exp -2Q 2 -2P 2 L t [AQ 2 + AP 2 



and 



Wna (q,p) 



T(n + a + 1/2) 



-Uq 2 /2 



:LT l/2 



- (q - i d ¥ ) 2 



L a-l/2 



2^ + i 



^ 2 + d^) 



-p 2 /2ZJ 



(19) 



respectively, where we have denoted irhWg (Q, P) and nhW na (q,p) as We(Q,P) and W na (q,p), 
respectively. 



3 Results and Discussion 

FIGs. 1 and 2 show, respectively, the WFs 1^02 and Wqs for the relative motion given by Eq. (fl9l as 
functions of dimensionless position q = q/l and momentum p = lp/h for two different dimensionless 
frequency values, U — u/ou, — 1 and 3. Contour plots showing the projections of the relevant WF 
onto (q, p) plane are also shown in these figures. In other words, each contour is a slice of given WF 
in the (q,p) plane. It should be noted that, while U — 1 corresponds to the case tJ = 0, u; = 3 
corresponds to switch on ZJ to the value oJq = 1, which causes localization in g[36l [37]. In addition 
to this pattern, derealization in p is observed in both figures. In other words, in FIG. 1(b) and FIG. 
2(b), there the dips and peaks of WFs in q are shifted towards the smaller q values compared with 
those in FIG. 1(a) and FIG. 2(a), whereas those of WFs in p are shifted towards higher p values. 

Having obtained a general expression for the WF of two interacting particles we now distinguish 
between the cases g = and g ^ 0, and deal with each case separately. This allows us to verify 
the consistency of the WFs obtained above with those found in the literature. In order to see this, 
we need to set g = first. In this case, we have only the solutions with (3 = +1/2 and —1/2 
corresponding to a = 1 and 0, respectively. We obtain, for n = 0, 1, 2, 3, . . ., 

7Tr , , v/?(2n+l)! f m 2 P 2 \ T ( 2 P 2 \ 

W nl (q, p) = , ; — ; exp — - — L 2n+1 [uJq 2 + — 20 
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which are the relative WFs corresponding to eigenvalues of the harmonic oscillators with 2n + 1 
eigenvalues. To obtain Eq. fl20l) . we have used the fact that every power of q 2 + d~ commutes with 
Laguerre polynomials with argument of q =)= % dp in Eq. ( [191) . and the identity 



- (q =f i dpf 



2 2n + 1 n\ 



H- 



2n+l 



(21) 



and we have derived a new identity between Hermite and Laguerre polynomials in the form of 

= (-l) n 2 n n\L n [2 (u 2 + v 2 )] e~ v \ (22) 



1 1 

H n (u+ -dv)H n (u- -&u)e 



The proof of Eq. fT2"2l can easily be done by using standard relations among these polynomials. In 
case of a = 0, we proceed as before by formally using Eqs. ([21]) and ([22]) to give 



W n0 (q, p) 



vr(2n)! 



n\ 2 2n T(n + 1/2) 



cxp 



ujq 



P 



1-^2 , V" 



L 2n [uq + _ 
2uj I \ uo 



(23) 



which are the relative WFs corresponding to eigenvalues of the harmonic oscillators with 2n eigen- 
values, again with n — 0, 1, 2, 3, A more general expression for these two cases can be found by 

noticing that arrangement of the coefficients in Eqs. (T20]) and ( 123]) leads to the pair of equations 



W n (q, p) 



exp 



9 — 9 

"~2 2UJ 



L 2n +i ( uq 2 + 



+ exp(-^--§.) I L ,,„(--'-i ^ 
or they may be combined into the form of 



W n (q,p) = (-l) n exp 



toq" jr \ T (—2,P 



2 ~z5 + *>' 



(24) 



The WFs constructed from the products of Eq. (T24l) with the Eq. (fT8l) define the WFs of two- 
noninteracting particles confined in a harmonic well potential in one dimension, or alternatively, 
they define WFs of a particle in a harmonic potential in two-space dimensions. When g ^ 0, which 
indicates that a would be greater than 1, then the total WF becomes 



W hn (q,p;Q,P) 
X U UQ 2 + 4P 



— 1)' exp 
■2\ n\ 



T (n + a + l/2) 



-2Q 2 - 2P 2 



-Uq 2 /2 
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(q - i d ¥ ) 



Lr l/2 



(q + * dp)" 



UJ 



-p 2 /2ZJ 



(25) 



Finally, by these considerations, we comment on W n (q, p) given by Eq. (T24|) . rather than Eq. 
(|25p . to better visualize the phase space behaviors of WFs presented in Figs. 1-2., i.e., how the 
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localization in q happens when the strength of spatial confinement u is increased. The use of the 
asymptotic expansion of the Laguerre polynomials for large order [38] yields Eq. (1241) to take form 



W n (q,p) 



n 



X COS 



(n + \)(ujq 2 + 1 ? = ) 

Z UJ 



-1/4 



1/2 



(26) 



from which one can easily find out where the WF is zero and where it takes negative values. For 
instance, Eq. (1261) has zeros when 



1^—2 . t 



2 uj 



1/2 



k 



0,±1,±2,..., 



(27) 



which clarifies the above comment that the localization in q appears as uj increases. Namely, the left 
hand side of Eq. (T2T1) can be rearranged, in dimensional units as usual, to give 



2H(q,p) 



7T 



UJ 



4(n +i) 



(k--) 2 h. 



(2f 



In fact, this last expression is called the symplectic area enclosed by an ellipse whose boundary is 
given by H(q,p) = E = p 2 /(2m) + uo 2 q 2 /2 and its minimum value is determined by the Gromov's 
non-squeezing theorem, i.e., 2TC(q,p)/uj > h [39]. Therefore, the projections of WFs onto (q,p) plane 



are elliptic energy shells whose eccentricity is given by e = sjl — (mw) 
the frequency uj = 1/m. 



and they are circles with 



In this paper, we introduced and solved the WFs of one-dimensional two particle Calogero- 
Sutherland system in which the particles obeying the Boltzman statistics interact mutually by the 
sum of quadratic and inversely quadratic pair potentials, and they are confined in an external har- 
monic potential as well. It is obvious that the technique introduced here can easily be extended to 
find explicit analytical expressions for WFs of 3-and N-body counterpart of the problem. Namely, 
by using Jacobi coordinates, after separating the center-of-mass coordinate, one can easily construct 
the remaining part of WF with N — 1 relative coordinates. Furthermore, due to the fact that, with 
particular choices of the coupling constant g, the radical equation of three dimensional isotropic 
oscillator and of hydrogen-like atom in both spherical and parabolic coordinates, one dimensional 
three body problem and the s-state of Morse potential[22] are all reduced to Calogero- Sutherland 
system, the results obtained here unify inherently the WFs of these quantum mechanical problems. 

As a final remark, we should point out that the attractive interaction, i.e., > g > —h 2 /8fi, is 
also present in the CSM. Therefore, one can easily compare the phase space behaviors of WFs of 
two different regimes as well. In particular, it should be noted that a particular choice of g, when 
> g > —h 2 /8fi, yields one dimensional band problem solved by Scarf [201 121] ■ This serves as a 
model for one dimensional dot arrays, as also indicated in the Introduction section. 
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(b) 

Figure 1: Three dimensional plots of WF W02 for the relative motion (Eq. ([19]) ) as a function 
dimensionless position q = q/l and momentum p = lp/h for dimensionless frequency (a) uj = uj/uj, 
1 and (b) uj = 3. 
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1 1 



W(b(«,p) 




(b) 

Figure 2: Three dimensional plots of WF W03 for the relative motion (Eq. (I19p ) as a function 
dimensionless position q = q/l and momentum p = Zp//i for dimensionless frequency (a) uj = uj/uj, 
1 and (b) cJ = 3. 
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